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Abstract: 

We search for CP violation in the decay — > KgK'^ using a data sample with 
an integrated luminosity of 977 fb~^ collected with the Belle detector at the KEKB e"^e~ 
asymmetric-energy collider. No CP violation has been observed and the CP asymmetry in 

— >■ KgK~^ decay is measured to be (—0.25 it 0.28 it 0.14)%, which is the most sensitive 
measurement to date. After subtracting CP violation due to — mixing, the CP 
asymmetry in D+ K^K+ decay is found to be (+0.08 ± 0.28 ± 0.14)%. 
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1 Introduction 

Studies of CP violation in charmed meson decays provide a promising opportunity to 
search for new physics beyond the standard model (SM) [1] in the absence of disagreement 
between experimental measurements and the SM interpretation of CP violation in K and 
B meson decays [2-4]. Recently, the LHCb collaboration has reported AAqp = (—0.82 it 
0.21 ± 0.11)% [5] where AAcp is the CP asymmetry difference between — t- K~^K~ ^ 
and — t- vr+vr" decays. Thereafter, the CDF collaboration has also announced /S.Acp = 
(— 0.62ib0.21 ±0.10)% [6], which strongly supports the non-zero AAcp measured from the 
LHCb collaboration. Together with results from the BaBar and Belle collaborations, the 
value of A. Ac p is significantly different from zero [7] . Taking into account that the indirect 
CP asymmetries in the two decays are approximately equal [8], AAqp can be expressed 
as 



where a^p and a™p denote direct and indirect CP violation, respectively, and {t)/T is 
the mean proper decay time of the selected signal sample in units of the lifetime [9]. 
The factor A{t)/T in eq. (1.1) depends on the experimental conditions and the largest 
value reported to date is 0.26 it 0.01 from the CDF measurement [6]. Therefore, AAcp 
reveals a significant direct CP violation difference between the two decays. Within the SM, 
direct CP violation in the charm sector is expected to be present only in singly Cabibbo- 
suppressed (SCS) decays, and even there is expected to be small, 0(0.1%) [10]. Hence, the 
current AAcp measurements engender questions of whether the origin of the asymmetry 
lies within [11-14] or beyond [15-18] the SM. The origin of AAcp calls for the precise 

^Throughout this paper, the charge-conjugate decay modes are imphed unless stated otherwise. 



AAcP = Aa^}. + a'SpMt)/r, 



(1.1) 
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measurements of Acp in — t- K^K^ and — )• vr^vr^. A complementary test is a precise 
measurement of Acp in another SCS charmed hadron decay, — t- as suggested 

in ref. [13]. As shown in figures 1(a) and 1(b), the decay — t- K^K~^ shares the same 
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Figure 1. Feynman diagrams of 13+ K^K^ and — > K^K decays. 

decay diagrams with — )• K^K~ by exchanging the spectator quarks, d ■(r^ u. Although 
there are additional contributions to the two decays as shown in figures 1(c) and 1(d), these 
are expected to be small due to helicity- and color-suppression considerations ^. Therefore, 
neglecting the latter contributions in Z)+ — )• K^K^ and — )• K^K~ decays, the direct 
CP asymmetries in the two decays are expected to be the same. 

In this paper, we report results from a search for CP violation in the decay — )• 
K^K^ that originates from — )• K^K^ decay, where decays to vr+vr^. The CP 
asymmetry in the decay, Aqp^ is then defined as 
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(1.2) 



where T is the partial decay width. In eq. (1.2), Aqp~^^°^^ is the CP asymmetry in the 
decay L>+ ^ K^K+ and A^p is that in —7- vr+TT decay induced by — mixing 
in the SM [19-21] in which the decay — t- vt+tt^ arises from — t- tt+tt^ together 



with a small contribution from — )• tt+tt , where the latter is known precisely from 
semileptonic decays, A^p = (-0.332 ± 0.006)% [2]. As shown in eq. (1.2), the product of 
the two small asymmetries is neglected. The decaying to the final state KgK~^ proceeds 
from D+ —7- K^K^ decay, which is SCS. In the SM, direct CP violation in SCS charmed 
meson decays is predicted to occur with a non-vanishing phase that enters the diagram 



^In helicity suppression, a spinless meson decaying to a back-to-back quark-antiquark pair is suppressed 
by the conservation of angular momentum. In color suppression, the final state quarks are required to carry 
the correct color charge in order for the final state to be colorless. 
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shown in figure 1(b) in the Kobayashi-Maskawa ansatz [22]. The current average of AAcp 
favors a negative value of direct CP violation in — t- K^K~ decay. Correspondingly, 
the CP asymmetry in — t- K^K^ decays is more likely to have a negative value since 
the two CP asymmetry terms shown in eq. (1.2) are negative. 



2 Methodology 

We determine A^p ^ by measuring the asymmetry in the signal yield 

-I'rec ~r JVrec 

where A'rec is the number of reconstructed decays. The asymmetry in eq. (2.1) includes 
the forward-backward asymmetry (Aps) due to 'y*-Z^ interference and higher order QED 
effects in e~^e~ — t- cc [23-25], and the detection efficiency asymmetry between and 
K~ {A^'^) as well as Acp- In addition, ref. [26] calculates another asymmetry source, 
denoted Ax>, due to the differences in interactions of and mesons with the material 
of the detector. Since we reconstruct the Kg with vr+vr" combinations, the tt~^tt~ detection 
asymmetry cancels out for Kg. The asymmetry of eq. (2.1) can be written as 

A^c^''^''^(cos0^r-,P^^+,cos0^^+,pfo) = A^s(cose^r-) 

+ Ar(pJ;i^^,cos0^b+) + ^i,(p^p (2.2) 

by neglecting the terms involving the product of asymmetries. In eq. (2.2), A^^p ^ is 
the sum of A^p~^^ and A^p as stated in eq. (1.2), where the former is independent of 
all kinematic variables while the latter is known to depend on the Kg decay time according 
to ref. [27], and Ap^ is an odd function of the cosine of the polar angle O'^f'^' of the 
momentum in the center-of-mass system (c.m.s.). A^^ depends on the transverse 
momentum Pp^+ and the polar angle of the K~^ in the laboratory frame (lab). Here, 
Ax> is a function of the lab momentum of the Kg. To correct for A^^ in eq. (2.2), 

we use the technique developed in our previous publication [28]. We use — )• K~tt~^ and 
Df —7- (pTr~^ decays where the (p is reconstructed with K^K~^ combinations and hence the 
K~^K~ detection asymmetry nearly cancels out [29] (the residual small effect is included 
in the systematic error). Since these are Cabibbo-favored decays for which the direct CP 
asymmetry is expected to be negligible, in analogy to eq. (2.2), A^^^~^^ and A^cc^'^^ 
can be written as 

4^ec"''""^(cos^^r-,p^l-,cos0^1,p^t+>cos^|,^) = A^l{cos O'^j^-'-) + Af" (p^^- , COS ; 



Afec""^" (cos0^r-,p;rt+,cos0j,^) = A^|j(cos0-r 



+ Af{p'^l^,cose':^), (2.3) 
A^|,(cos0-?-) 

Af{p'^l^,cose':^). (2.4) 
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Thus, with the additional term in A^^~^^ one can measure by subtracting 
A^c^'^^ from Aj!^^~^^ assuming the same Aps for D*^ and Df mesons. We also 
obtain Ax> according to ref. [26]. After these A^^ and Ax> corrections ^, we obtain 

^,^,^^^5^— (cos0^?-^-) = Agp^"^^"^^ + ^^s(cose^?-^-)- (2.5) 

We subsequently extract Aqp and App as a function of cos0^™'^' by taking sums and 
differences: 

,D+^KOK+. A^c (+ COS g^T-") + ^gc (- cos O'j^T) 

A^p 5 (Icos^^'f'^i) = ^ 2 , 

(2.6a) 

^^^(1 cos^-ri) = — ^+ \ ^ \ (2.6b) 

Note that extracting Acp in eq. (2.5) using eq. (2.6a) is crucial here to cancel out the Belle 
detector's asymmetric acceptance in cos^^™''^'. 

3 Data and event selections 

The data used in this analysis were recorded at the T{nS) resonances (n = 1,2,3,4,5) or 
near the T(4S) resonance with the Belle detector at the e~^e~ asymmetric-energy collider 
KEKB [30]. The data sample corresponds to an integrated luminosity of 977 fb^^. The 
Belle detector is a large solid angle magnetic spectrometer that consists of a silicon vertex 
detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprising CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return 
located outside the coil is instrumented to detect mesons and to identify muons. A 
detailed description of the Belle detector can be found in ref. [31]. 

Except for the tracks from Kg decays we require charged tracks to originate from the 
vicinity of the interaction point (IP) by limiting the impact parameters along the beam 
direction (z-axis) and perpendicular to it to less than 4 cm and 2 cm, respectively. All 
charged tracks other than those from Kg decays are identified as pions or kaons by requiring 
the ratio of particle identification likelihoods, Ck /{^k+^-k), constructed using information 
from the CDC, TOF, and ACC, to be larger or smaller than 0.6, respectively [32]. For both 
kaons and pions, the efficiencies and misidentification probabilities are about 90% and 5%, 
respectively. 

We form Kg candidates adopting the standard Belle Kg criteria [33], requiring the 
invariant mass of the charged track pair to be within [0.4826, 0.5126] GeV/c^. The "loose" 
Kg candidates not satisfying these standard selections are also used in this analysis with 
additional requirements described later. 

The Kg and K~^ candidates are combined to form a candidate by fitting their 
tracks to a common vertex; the candidate is fitted to the independently measured IP 



^We define A^'^ = [N''^ - N*" ]/[N''^ + iV'^ ]. Hence A*" = -A^* . 
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profile to give the production vertex. To remove combinatorial background as well as 
mesons that are produced in possibly CP- violating B meson decays, we require the 
meson momentum calculated in the c.m.s. to be greater than 2.5 and 3.0 GeV/c for 

the data taken at the T(4S') and T{5S) resonances, respectively. For the data taken below 
T(4S'), where no B mesons are produced, we apply the requirement p^+>2.0 GeV/c. In 
addition to the selections described above, we further optimize the signal sensitivity with 
four variables: the goodness-of-fit values of the decay- and production- vertex fits Xd 
and Xp7 the transverse momentum of the in the lab Pj^^+, and the angle ^ between the 

momentum vector (as reconstructed from its daughters) and the vector joining the 
production and decay vertices. We optimize the requirement on these four variables with 
the standard and loose Kg selections by maximizing Ms / V-^S + -N^B , where Ms + Mb and 
Mb are the yields in the K^K~^ invariant mass signal ([1.860, 1.884] GeV/c^) and sideband 
([1.843,1.855] and [1.889,1.901] GeV/c^) regions, respectively. The optimal set of (x|,, 
Xp, Vt\+7 6 requirements are found to be (<100, <10, >0.30 GeV/c, <40°), (<100, 
<10, >0.25 GeV/c, <115°), and (<100, <10, >0.20 GeV/c, <125°) for the data taken 
below the T(4S'), at the T(45'), and at the T(55), respectively. Note that is highly 
correlated with and ^; hence, a tighter requirement on the T(5S') sample results 

in looser p^^+ and ^ requirements and vice versa for the data taken below the T(45). The 

candidates with the loose requirement are further optimized with two additional 
variables: the of the fit of tracks from the decay and the kaon from the meson 
decay to a single vertex {Xxkh) ^^'^ the angle C between the Kg momentum vector (as 
reconstructed from its daughters) and the vector joining the and Kg decay vertices. The 
two variables are again varied simultaneously and the optimum is found to be Xkhh^^ 
C<3° for all data. The inclusion of candidates with the loose Kg requirement improves 
the statistical sensitivity by approximately 5%. After the final selections described above, 
we find no significant peaking backgrounds — for example, -D"*" — t- Tr~^TT~K^ decays — in the 
Monte Carlo (MC) simulated events [34]. Figure 2 shows the distributions of M[KgK^) 
and M{KgK^) together with the results of the fits described below. 

Each — )• KgK^ signal is parameterized as two Gaussian distributions with a 
common mean. The combinatorial background is parameterized with the unnormalized 
form eO^+I^Mi^s^"^) ^ where a and /? are fit parameters. The asymmetry and the sum of 
the D'^ and yields are directly obtained from a simultaneous fit to the and D" 
candidate distributions. Besides the asymmetry and the sum of the D'^ and D~ yields, 
the common parameters in the simultaneous fit are the widths of the two Gaussians and 
the ratio of their amplitudes. The asymmetry and the sum of the — )• KgK~^ and 
D- K^K- yields from the fit are (+0.048 ± 0.275)% and 276812 ± 1156, respectively, 
where the errors are statistical. 

In order to measure the CP asymmetry in — )■ KgK~^ decays, we must also re- 
construct — K^-K^ and Df — (pir^ decays: see eqs. (2.2), (2.3), and (2.4). For 
the reconstruction of the — t- K^tt^ and D'^ — t- (pir'^ decays, we require the same track 
quality, particle identification, vertex fit quality, and requirements as used for the recon- 
struction of the — )• KgK~^ decays, where the mass window for the (p is ±16 MeV/c^ [29] 
of the nominal (p mass [2]. 
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Figure 2. Distributions of M{Kl.K+) (left) and M{K'^K-) (right). Dots are the data while the 
histograms show the results of the parameterizations of the data. Open histograms represent the 
— >■ KgK"^ signal and shaded regions are combinatorial background. 
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Extraction of Acp in the decay — )► KgK^ 



To obtain Af-^ , we first extract A^c '^'^'^ from a simultaneous fit to the mass distributions 
of Df and D~ candidates with similar parameterizations as for — )• K^K^ decays 
except that, for the Df — )• (pn^ signal description, a single Gaussian is used. The values 
of A^c^'^^ are evaluated in 10x10x10 bins of the three-dimensional (3D) phase space 
cos 9]^+, cos0^+-^-). Each — ^ K^-ir^ and Z)° — ^ K^ir" candidate is then weighted 

with a factor of 1 — A^c ^'^'^ and 1 + A^cc ^"^'^ , respectively, in the corresponding bin of 
this space. After this weighting, the asymmetry in the — )■ K~i\:'^ decay sample becomes 
A^ . The detector asymmetry, Af , is measured from simultaneous fits to the weighted 
M{K^'K^) distributions in 10x10 bins of the 2D phase space (p)^^_, cos6^^_) with similar 
parameterizations as used for — )■ K^K~^ decays except that, for the — )• K~t:~^ 
signal description, a sum of a Gaussian and bifurcated Gaussian is used. Figure 3 shows 
the measured A^ in bins of and cos 6^^_ together with A^^^~^^ for comparison; 

we observe that A^^,"^^ '^^ shows a cos 9^^_ dependency that is inherited from A^^ while 
Af~ does not. The average of Af over the phase space is (—0.150 it 0.029)%, where the 
error is due to the limited statistics of the — )• K~'a~^ sample. 

Based on a recent study of A-p [26], we obtain the dilution asymmetry in bins of 
lab momentum. For the present analysis, Ad is approximately 0.1% after integrating over 
the phase space of the two-body decay. 

The data samples shown in figure 2 are divided into 10x10x16 bins of the 3D phase 
space {p^^\+ , cos , ) . Each — )• K^K^ candidate is then weighted with a factor 

of (1 ^ ^)(1 T Av) in this space. The weighted M{K^K^) distributions in bins of 
cosO^fy^'^' are fitted simultaneously to obtain the corrected asymmetry. We fit the linear 
component in cos^^™'*^' to determine Afb\ the Acp component is uniform in cos^^™'*^'. 

Figure 4 shows A^p^^^^ and A^^ as a function of | cos 9^'^' \ . From a weighted average 
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Figure 3. The map in bins of pJf and cosS'**" of the K obtained with the ^ K n 



and _D+ 07r+ samples (triangles). The A 



map is also shown (rectangles). 



over the | cos 9'^'^' \ bins, we obtain A^p « = (-0.246 ± 0.275)%, where the error is 
statistical. 

5 Systematic uncertainty 

The entire analysis procedure is validated with fully simulated MC events [34] and the 
result is consistent with null input asymmetry. We also consider other sources of systematic 
uncertainty. The dominant one in the Acp measurement is the A^^ determination, the 
uncertainty of which is mainly due to the statistical uncertainties in the — )• K~Tr~^ 
and — )• (j)7r~^ samples. These are found to be 0.029% and 0.119%, respectively, from a 
simplified simulation study. A possible Acp in the — t- K'tt^ final state is estimated 
using Acp = —y sin 6 sin (f)^/R [35]. A calculation with 95% upper and lower limits on 
mixing and CP violation parameters y, (f>, and strong phase difference 5 and 
Cabibbo suppression factor R from ref. [3], Acp in the — )• K~tt~^ final state is estimated 
to be less than 0.005% and this is included as one of systematic uncertainties in the A^^ 
determination. As reported in our previous publication [29], the magnitude of A^^ for 
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Figure 4. Measured Aqp (top) and Aps (bottom) values as a function of | cos6'^f '^ l. In the top 
plot, the dashed line is the mean value of Acp while the hatched band is the ildtotai interval, 
where Ctotai is the total uncertainty. 



the cj} reconstruction in Df — )• (j)TT~^ decays is 0.051%, which is also added to the systematic 
uncertainty in the A^^ measurement. By adding the contributions in quadrature, the 
systematic uncertainty in the A^^ determination is estimated to be 0.133%. We estimate 
0.008% and 0.021% systematic uncertainties due to the choice of the fitting method and 
that of the cosO'j^'^' binning, respectively. Finally, we add the systematic uncertainty in 
the Ajy correction, which is 0.010% based on ref. [26]. The quadratic sum of the above 
uncertainties, 0.135%, is taken as the total systematic uncertainty. 



6 Results 

We find A^p ^ = (—0.246 ± 0.275 ± 0.135)%. This measurement supersedes our 

previous determination of A^p ^ [28]. In Table 1, we compare all the available mea- 
surements and give their weighted average. 

According to Grossman and Nir [27], we can estimate the experimentally measured 
CP asymmetry induced by SM — mixing, Aqp. The efficiency as a function of Kg 
decay time in our detector is obtained from MC simulated events. The efficiency is then 
used in eq. (2.10) of ref. [27] to obtain the correction factor that takes into account, for 
A^p, the dependence on the kaon decay time. The result is 0.987 ± 0.007. By multiplying 
the correction factor 0.987 it 0.007 and the asymmetry due to the neutral kaons [2], we find 
the experimentally measured A^p to be (—0.328 it 0.006)%. 
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Experiment 


Acp ' (%) 


FOCUS [36] 
CLEO [37] 

Belle (this measurement) 


+7.1 ±6.1 ±1.2 
-0.2 ± 1.5 ±0.9 
-0.246 ±0.275 ±0.135 


New world average 


-0.23 ±0.30 



Table 1. Summary of A^p ^ measurements (where the first uncertainties are statistical and 
the second systematic), together with their average (assuming the uncertainties to be uncorrelated, 
the error on the average represents the total uncertainty). 

7 Conclusion 

We report the most sensitive CP asymmetry measurement to date for the decay — )• 
K^K~^ using a data sample corresponding to an integrated luminosity of 977 fb~^ collected 
with the Belle detector. The CP asymmetry in the decay is measured to be (—0.25 ± 
0.28 ± 0.14)%. After subtracting the contribution due to - mixing {A^p), the CP 
asymmetry in the charm decay {A^^~^^^^*) is measured to be (±0.08 ± 0.28 ± 0.14)%, 
which can be compared with direct CP violation in — t- K'^K~ . For the latter the 
current averages of l^Acp and CP asymmetry in — )• K~ favor a negative value [3]. 
Our result, on the other hand, does not show this tendency for — >■ K^K~^ decays, albeit 
with a significant statistical uncertainty. 
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